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METHOD AND APPARATUS FOR CONTROLLING 



A PLURALITY OF LOCOMOTIVES 



BACKGROUND 

[0001] The present invention relates generally to the field of controlling 
locomotives and more specifically to the field of controlling a plurality of locomotives 
so as to optimize an objective function while satisfying constraints on a set of 
performance parameters. 

[0002] Freight trains are often hauled by multiple locomotive ensembles 
("consists") placed together at the front or rear of the train or dispersed among the 
freight cars. A single crew at the front of the train coordinates all the locomotive 
throttles and brake commands via a cable called the multiple unit line ("MU-line") 
that runs among the locomotives. The MU-line can also be implemented virtually 
using radio links. 

[0003] For historical reasons, locomotive throttles are typically set up to have 
discrete settings ("notches," typically designated "idle" and "Nl" through "N8") 
corresponding to different levels of pulling power. As locomotive speed varies, on 
board electronic control circuitry regulates the power production by electric traction 
motors to be approximately constant at a corresponding level defined for each notch. 
The locomotive carrying the train driver is designated the "lead locomotive"; the other 
locomotives are said to be "trailing" or "in trail." As conventionally coordinated, 
each locomotive in trail is operated at the same notch set by the train driver for the 
lead locomotive. As the train driver manipulates the throttle in the lead locomotive, 
the trailing locomotives conventionally follow in lock step producing a range of 
horse-power from idle up to the combined rating of the units, track conditions and 
speed limits permitting. 

[0004] For each individual locomotive, fuel efficiency, as measured for example in 
pounds of fuel per horsepower-hour, generally improves with increasing power 
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output. When less than full power is required from the overall consist, therefore, 
there is an efficiency penalty that results from operating all locomotives at the same 
low power setting. For example, a typical three locomotive consist operating with all 
locomotives in N4 produces substantially the same power but consumes more fuel 
than the same consist operating with one locomotive in N8 and the other two 
locomotives in idle. Opportunities exist, therefore, to improve overall consist 
efficiency by allowing the notches to be set independently on each locomotive. 

SUMMARY 

[0005] The opportunities described above are addressed, in one embodiment of the 
present invention, by an apparatus for controlling a plurality of locomotives, the 
locomotives being responsive to a plurality of discrete actual commands, the 
apparatus comprising: a combination generator adapted for generating combinations 
of the discrete actual commands to yield a command combination set; a performance 
calculator adapted for calculating a performance parameter set from the command 
combination set; a feasible combination selector adapted for selecting a feasible 
combination subset of the command combination set as a function of a discrete 
performance setpoint, a performance tolerance, and the performance parameter set; an 
objective function calculator adapted for calculating an objective function set from the 
feasible combination subset; and an optimal command selector adapted for selecting 
an optimal command combination from the feasible combination subset 
corresponding to an optimum value of the objective function set. 

[0006] The present invention is also embodied as a method for controlling a 
plurality of locomotives, the locomotives being responsive to a plurality of discrete 
actual commands, the method comprising: generating combinations of the discrete 
actual commands to yield a command combination set; calculating a performance 
parameter set from the command combination set; selecting a feasible combination 
subset of the command combination set as a function of a discrete performance 
setpoint, a performance tolerance, and the performance parameter set; calculating an 
objective function set from the feasible combination subset; and selecting an optimal 
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command combination from the feasible combination subset corresponding to an 
optimum value of the objective function set. 

DRAWINGS 

[0007] These and other features, aspects, and advantages of the present invention 
will become better understood when the following detailed description is read with 
reference to the accompanying drawings in which like characters represent like parts 
throughout the drawings, wherein: 

[0008] Figure 1 illustrates a block diagram in accordance with one embodiment of 
the present invention. 

[0009] Figure 2 illustrates a block diagram in accordance with a more specific 
embodiment of the embodiment of Figure 1. 

[0010] Figure 3 illustrates a block diagram in accordance with a more specific 
embodiment of the embodiment of Figure 2. 

[0011] Figure 4 illustrates a block diagram in accordance with another more 
specific embodiment of the embodiment of Figure 2. 

[0012] Figure 5 illustrates a block diagram in accordance with still another more 
specific embodiment of the embodiment of Figure 2. 

[0013] Figure 6 illustrates a block diagram in accordance with another 
embodiment of the present invention. 

DETAILED DESCRIPTION 

[0014] In accordance with one embodiment of the present invention, Figure 1 
illustrates a block diagram of an apparatus 100 for controlling a plurality of 
locomotives responsive to a plurality of discrete actual commands. Apparatus 100 
comprises a combination generator 110, a performance calculator 130, a feasible 
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combination selector 150, an objective function calculator 170, and an optimal 
command selector 190. In operation, combination generator 110 generates 
combinations of the discrete actual commands to yield a command combination set 
120. From command combination set 120, performance calculator 130 calculates a 
performance parameter set 140. Feasible combination selector 150 then selects a 
feasible combination subset 160 from command combination set 120 as a function of 
a discrete performance setpoint 152, a performance tolerance 154, and performance 
parameter set 140. From feasible combination subset 160, objective function 
calculator 170 calculates an objective function set 180. Optimal command selector 
190 then selects an optimal command combination 200 from feasible combination 
subset 160 corresponding to an optimum value of objective function set 180. 

[0015] In a more detailed embodiment in accordance with the embodiment of 
Figure 1, combination generator 110 generates all possible combinations of discrete 
actual commands. Such an embodiment provides a "full optimization" by exhaustion. 

[0016] In an alternative more detailed embodiment in accordance with the 
embodiment of Figure 1, combination generator 110 generates only combinations 
wherein discrete actual commands corresponding to selected locomotives are equal. 
In other words, the selected locomotives are operated in lock step while the remaining 
locomotives are free to be operated independently. Such an embodiment provides a 
"partial optimization." For some applications, the partial optimization is significantly 
less expensive to implement than the full optimization and is therefore preferable 
despite the sacrifice of some optimization benefit. 

[0017] In another more detailed embodiment in accordance with the embodiment 
of Figure 1, performance calculator 130 calculates total power of the plurality of 
locomotives. In such an embodiment, discrete performance setpoint 152 is a request 
for a level of total consist power, performance tolerance 154 is a power tolerance, and 
the discrete actual commands are typically notches. Feasible combination subset 160 
therefore contains only notch combinations that result in a total consist power within 
the power tolerance of the request. 
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[0018] In still another more detailed embodiment in accordance with the 
embodiment of Figure 1, discrete performance setpoint 152 has a finer resolution than 
the discrete actual commands. As used herein, a discrete signal X "has a finer 
resolution" than a discrete signal Y if the number of discrete values of discrete signal 
X is greater than the number of discrete values of signal Y. Consider, for example, a 
consist of three locomotives each responsive to nine levels (i.e., idle plus Nl through 
N8). If operated conventionally in lock step, only nine discrete values of total consist 
power can be achieved. If operated in accordance with the present embodiment, 
however, at least twenty-five discrete values of total consist power can be achieved. 

[0019] In yet another more detailed embodiment in accordance with the 
embodiment of Figure 1, objective function calculator 170 calculates a quantity 
selected from the group consisting of fuel efficiency, fuel remaining in each 
locomotive, engine wear, cab noise, and deviation from discrete performance setpoint 
152. 

[0020] In accordance with a more specific embodiment of the embodiment of 
Figure 1, Figure 2 illustrates a block diagram wherein performance calculator 130 
further comprises an observer 210 and a performance module 212. In operation, 
observer 210 estimates a set of estimated state variables 220 corresponding to the 
locomotives, and performance module 212 calculates performance parameter set 140 
from command combination set 120 and estimated state variables 220. 

[0021] Observer 210 comprises an update module 214 and a process model 252. 
In operation, update module 214 generates state update signals 230 from measured 
process outputs 240 and predicted process outputs 250. Process model 252 
propagates estimated state variables 220 and calculates predicted process outputs 250 
using state update signals 230 and a set of measured process inputs 260. Examples of 
observer 210 include, without limitation, Luenberger observers, extended observers, 
and finite- or infinite-horizon Kalman filters, and extended Kalman filters. 
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[0022] In a more detailed embodiment in accordance with the embodiment of 
Figure 2, observer 210 is an extended Kalman filter. 

[0023] In accordance with a more specific embodiment of the embodiment of 
Figure 2, Figure 3 illustrates a block diagram wherein process model 252 comprises a 
cylinder charge estimator 270, a cylinder frequency calculator 300, a multiplier 330, 
an electrical power calculator 350, and an integrator 370. In operation, cylinder 
charge estimator 270 estimates a cylinder charge 290 (in units, for example, of liters 
per deflagration) from a fuel command 280 using a cylinder charge versus fuel 
command table 272. Cylinder frequency calculator 300 calculates a cylinder 
frequency 310 (in units, for example, of deflagrations per second) from an engine 
speed measurement 320. To yield a fuel flow estimate 340 (in units, for example, of 
liters per second), multiplier 330 multiplies cylinder charge 290 by cylinder frequency 
310. Electrical power calculator 350 then calculates an electrical power estimate 360 
from fuel flow estimate 340 using an electrical power versus fuel flow table 352. Fuel 
flow estimate 340 is integrated with respect to time by integrator 370 to yield a 
consumed fuel volume estimate 380. In the Figure 3 embodiment of process model 
252: state variables 220 comprise elements of cylinder charge versus fuel command 
table 272, elements of electrical power versus fuel flow table 352, and consumed fuel 
volume estimate 380; measured process inputs 260 comprise fuel command 280 and 
engine speed measurement 320; predicted process outputs 250 comprise electrical 
power estimate 360 and consumed fuel volume estimate 380; and measured process 
outputs 240 comprise an electrical power measurement and a consumed fuel volume 
measurement (not shown). 

[0024] In accordance with another more specific embodiment of the embodiment 
of Figure 2, Figure 4 illustrates a block diagram wherein process model 252 
comprises a fuel flow calibration module 420, an electrical power calculator 350, and 
an integrator 370. In operation, fuel flow calibration module 420 generates a fuel 
flow estimate 340 from a fuel flow measurement 410 using a fuel flow calibration 
table 422. From fuel flow estimate 340, electrical power calculator 350 calculates an 
electrical power estimate 360 using an electrical power versus fuel flow table 352. 
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Integrator 370 then integrates fuel flow estimate 340 over time to yield a consumed 
fuel volume estimate 380. In the Figure 4 embodiment of process model 252: state 
variables 220 comprise elements of fuel flow calibration table 422, elements of 
electrical power versus fuel flow table 352, and consumed fuel volume estimate 380; 
measured process inputs 260 comprise fuel flow measurement 410; predicted process 
outputs 250 comprise electrical power estimate 360 and consumed fuel volume 
estimate 380; and measured process outputs 240 comprise an electrical power 
measurement 390 and a consumed fuel volume measurement (not shown). 

[0025] In accordance with still another more specific embodiment of the 
embodiment of Figure 2, Figure 5 illustrates a block diagram wherein process model 
252 comprises a fuel flow calculator 430 and an integrator 370. In operation, fuel 
flow calculator 430 generates a fuel flow estimate 340 from an electrical power 
measurement 390 using an electrical power versus fuel flow table 352. Integrator 370 
integrates fuel flow estimate 340 over time to yield a consumed fuel volume estimate 
380. In the Figure 5 embodiment of process model 252, state variables 220 comprise 
elements of electrical power versus fuel flow table 352 and consumed fuel volume 
estimate 380; measured process inputs 260 comprise electrical power measurement 
390; predicted process outputs 250 comprise consumed fuel volume estimate 380; and 
measured process outputs 240 comprise a consumed fuel volume measurement (not 
shown). 

[0026] In accordance with another embodiment of the present invention, Figure 6 
illustrates a block diagram wherein apparatus 100 further comprises a transition 
control module 440. In operation, transition control module 440 computes an applied 
command combination 450 from discrete performance setpoint 152 and optimal 
command combination 200 using a transition control method. The transition control 
method comprises ramping elements of applied command combination 450 toward 
respective elements of optimal command combination 200, but delaying the ramping 
of a subset of the elements to minimize an adverse performance response. As used 
herein, an "adverse performance response" refers to an undesirable characteristic of a 
time history of a performance parameter. For example, a response to a requested 
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increase in total consist power wherein the time history of total consist power was not 
monotonically increasing would be an adverse performance response. 

[0027] By way of example, but not limitation, suppose: discrete performance 
setpoint 152 is a request for a level of total consist power; the discrete actual 
commands are notches; there are two locomotives in the consist; the leading 
locomotive is currently at N8 and the trailing locomotive is currently at idle, 
represented as the ordered notch pair (8, 0); and the next value of discrete 
performance setpoint 152 corresponds to notch pair (5, 5). If the notches for both 
locomotives are ramped together, a first sequence would be (7, 1), (6, 2), (5, 3), (5, 4), 
(5, 5). If the ramp of the lead locomotive is delayed two periods, a second sequence 
would be (8, 1), (8, 2), (7, 3), (6, 4), (5, 5). For locomotives with notches equally 
spaced in power, the first sequence exhibits an adverse performance response in that 
the total consist power fails to increase in the first two time periods. The second 
sequence minimizes that adverse performance response by causing the total consist 
power to increase immediately. More complicated situations arise when the notch 
powers are not equally spaced. 

[0028] In a more specific embodiment in accordance with the embodiment of 
Figure 6, the act of delaying ramping further comprises waiting a specified time. The 
specified time is a function of optimal command combination 200 and of a previous 
value of applied command combination 450. In some embodiments, the specified 
time is a tabulated function pre-computed to minimize an adverse performance 
response. 

[0029] In another more specific embodiment in accordance with the embodiment 
of Figure 6, the act of delaying ramping further comprises calculating a predicted 
performance response due to ramping elements of applied command combination 
450, and delaying ramping of a subset of the elements until the predicted performance 
response satisfies a transition performance criterion. In contrast to the previous 
embodiments, instead of waiting a pre-computed time, this embodiment decides in 
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real time whether the response to ramping applied command combination 450 will be 
acceptable. 

[0030] All of the above described elements of embodiments of the present 
invention may be implemented, by way of example, but not limitation, using singly or 
in combination any electric or electronic devices capable of performing the indicated 
functions. Examples of such devices include, without limitation: analog devices; 
analog computation modules; digital devices including, without limitation, small-, 
medium-, and large-scale integrated circuits, application specific integrated circuits 
(ASICs), and programmable logic arrays (PLAs); and digital computation modules 
including, without limitation, microcomputers, microprocessors, microcontrollers, and 
programmable logic controllers (PLCs). 

[0031] In some implementations, the above described elements of the present 
invention are implemented as software components in a general purpose computer. 
Such software implementations produce a technical effect of controlling the plurality 
of locomotives so as to optimize a selected objective function. 

[0032] While only certain features of the invention have been illustrated and 
described herein, many modifications and changes will occur to those skilled in the 
art. It is, therefore, to be understood that the appended claims are intended to cover 
all such modifications and changes as fall within the true spirit of the invention. 
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